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Abstract 

A method of solving the eikonal equation, in either flat or curved 
space-times, with arbitrary Cauchy data, is extended to the case of 
data given on a characteristic surface. We find a beautiful relationship 
between the Cauchy and characteristic data for the same solution, 
namely they are related by a Legendre transformation. From the 
resulting solutions, we study and describe the wave-front singularities 
that are associated with their level surfaces (the characteristic surfaces 
or "big wave- fronts" ) . 



1 Introduction 



The high frequency limit of the wave equation is given by the eikonal equa- 
tion, written in an arbitrary space-time as: 

g ab {x a )d a S{x a )d b S{x a ) = (1) 

where the x a = (x l ,t) are any local coordinates, and g ab ( x a ) is the met- 
ric of the given space-time. The level surfaces of a solution of Eq. ([!]), 
S(x a ) = const., (which need not be smooth every place and could have self- 
intersections) are three dimensional characteristic surfaces (the "big wave- 
fronts" in the terminology of Arnold |L|), and the sections t = constant of 
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these surfaces are the two-dimensional ("small") wave fronts. The vector 
field l a = g ab dbS is tangent to the null geodesic that generate the character- 
istic surfaces. 

In flat space-time the eikonal equation becomes 

V ab d a Sd b S = (d t S) 2 - (d x S) 2 - (d y S) 2 - (d z S) 2 = 0. (2) 

In Sec. 2 we review the method @ to give a general solution of the eikonal 
equation in flat space-time adapted to appropriate Cauchy data given at t — 
t . In Sec. 3 we modify the method so that the eikonal equation is solved with 
arbitrary characteristic data given at null infinity. The relation between both 
methods is studied in Sec. 4 where we find that the Cauchy and characteristic 
data are related by a Legendre transformation. The wave-front singularities 
of the level surfaces of the resulting solutions are described parametrically 
in Sec. 5 and finally the generalization of our results to asymptotically flat 
spaces-times is given in Sec. 6. 

2 Solutions of the Eikonal Equation 

^From the point of view of the theory of partial differential equation the 
eikonal equation is a homogeneous first order non-linear partial differential 
equation; there exist a solution S*, called 'complete integral', depending on 
three arbitrary const ant s||, e.g., in flat space-time the function 

S*(x i ,t,a i )=x i a i -t^2(a i ) 2 (3) 

is easily seen to satisfy (0). 

Remark 1. The fact that the equation is homogeneous plays no role in this 
section, but will be crucial in the generalization to characteristic data. 

Remark 2. From now on we will treat the problem of the eikonal in flat 
space; we leave the discussion of general space-times to the end. 

It is possible to generate a "general integral", i.e., a solution of the eikonal 
equation depending on an arbitrary function, by means of the following pro- 
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cedure: First, define the function S'**(x J , t, aj) of the coordinates and the free 
parameters oti as 



S**(x\ t, a t ) = S*(x\ t, a,) - H{ ai ), (4) 

where H(a.i) is any function of the a^s. 

Next, think of the a^s as functions of the space-time points obtained from 
the following conditions: 

dS**/d(Xi = dS*/don - dH/dai = (5) 

Equation (^) determine the space-time dependence of the a^s (c^ = Ai(x % , t)) 
given that 

d 2 S**(x\t, ai ) 



daidctj 



+ 0. (6) 



This condition can fail in lower dimensional regions called the caustics. This 
issue will be returned to in Sec. 5. 

Substituting cti = Ai(x l ,t) into equation (f|) we eliminate the a i: and 
obtain 

S**(x\ t) = S*(x\ t, A(x\ t)) - H(A t (x\ £)). (7) 
It is easy to verify that, because of the condition (|]), 

d a S** = d a S* (8) 

which means that S**(x l ,t) is a new solution of the eikonal equation (TO) 
determined by an arbitrary function H. We can determine the free function 
H so that the solution (|7|) satisfies initial Cauchy data at t = t . We denote 
the Cauchy data by Sc a uchy(% 1 )- Conditions (||) imply that at t = t , oti = 
dScauchyl dx l . Inverting these relations we obtain x % = and replacing 

them in (0) at t = t we find the sought for relation: 

H(a t ) = S*{X\ ai )^ ai ) - S Cauchy (X\a t )). (9) 

This last equation relates the arbitrary function H(ai) with the Cauchy 
data, Scauch y (x*), at t = to an d allows us to construct solutions of the eikonal 
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equation in flat space-time for any initial data. 

We now change the set of the a l, s for new parameters that are more appro- 
priate to the study of asymptotically flat spaces and we rewrite our previous 
equations in terms of them. A complete integral of the eikonal equation, 
equation (1), can be written in terms of new parameters (/3, (, () as: 

s*(x a ,p,c,C) = Px a Uc,& (io) 

where 

UC, = * -=r ((i + CO, (C + 0, -<(C - 0, (CC - 1)) (ii) 
v2(l + CC) 

is the null covector pointing in the (C,C) direction. The (C? C) are the stere- 
ographic coordinates that parametrize the sphere of null directions. 

/,From ([IT]) and (|H]), we get the relations between the new parameters 
(P, C, C), and the old a's: 



a i 



«2 



P c + c 

^i + CC 

._g c-c 
l V2i + a 



a3 = TITTcc' (12) 



and 



In terms of the new parameters Eq. (ffl) reads: 

S**(x a , (3, c, C) = /^X(C, - H(/3, C, 0, (14) 

where H(/3,£, C) is an arbitrary function that will be determined by the initial 
conditions. Conditions © on (/5, C ,C) become 



fjTT 

x a UU)--Qp(P,C,()\a = o (is) 
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(3x a m a ((,0- m(P,(,()y = (16) 
[3x a m a ((X)-m(f3,(,0\p = 0. (17) 

Remark 3. We have replaced the derivatives with respect to ( and ( respec- 
tively by 

3 = (! + CC)| 

and used the fact that S/ a (C C) = m a and 9 9Z a = n a —l a where (l a , n a , m a , m a ) 
form a null Minkowski space tetrad for each ((,()■ 

The function H(j3,(, C) can De determined by means of the same proce- 
dure using the conditions c*i = dScauchy/dx 1 at t — to and relations ([12]) to 
obtain the x % = X l (/3, (, (), and finally rewriting (|) 

H((3, C, = S*(X\(3, C, 0^0, P, C, C) - ScauchyiX'iP, C, C))- 

3 Characteristic Data for the Eikonal Equa- 
tion 

The eikonal equation, being hyperbolic, admits a characteristic formulation. 
Even though the results of this section can be applied to any characteristic 
hypersurface in Minkowski, in flat space (as in any asymptotically simple 
space-time) there are two preferred characteristic surfaces, namely future 
and past null infinity respectively. In the following we will formulate the 
characteristic problem in terms of data given at future null infinity, 3 + . 3 + 
has the topology of S 2 x R ; we choose Bondi coordinates on it, namely (£, 
() on the S 2 and the retarded time ub along R. In an analogous manner as 
for the Cauchy problem, the characteristic data at future null infinity will be 
defined by a function of C> C )j 

S Characteristic = L(ub, C> C)- (18) 
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The goal of this section is to develop a method to construct solutions of 
the eikonal equation geometrically adapted to the characteristic data, Eq.( 
0), at 3 + . 

Remark 4. In asymptotically flat space-times, in the neighborhood of future 
null infinity, J + , there is a preferred class of coordinates referred to as Bondi 
coordinates. Given a Bondi system (ub,CX) a ^ 3 + , a new system^ (u,(,() 
can be defined by (u,(X) = (L(ub, (, C)> d 0- The characteristic data, 
can be thought of as being generated by this coordinate change representing a 
one parameter family of arbitrary u = const, slices of3 + . 

In flat space-times we can define a two parameters family of null surfaces 

by: 

S* = x a l a ((,() (19) 

As was pointed out in Eq.(|19|) has a dual interpretation. For (C,C) 
kept constant, its level surfaces define null planes intersecting the time axis 
at a time equal the value of S* and with its direction given by ((, Q; on the 
other hand, for a fixed value of x a it represents the light cone cut at 3 + of 
the space time point x a in the interior, i.e., it represents the intersection of 
the null cone from x a with 3 + . 

We can think of the characteristic data ( |T8D geometrically, as defining a 
one parameter family of cuts at 3 + in terms of u = L(ub, C> C) = const. (It is 
assumed that this can be inverted so that the cuts are given by u B = L~ l {u = 
const, (,()■) With this point of view, we construct a solution of the eikonal 
equation (|l|) (corresponding to the characteristic data L(ub,(X)), such that 
the family of null surfaces in the interior are defined by the null geodesies 
normal to the family of cuts at infinity given by L(ub,(>0 = const.. In 
order to do so we will generalize the method of section II. 

Defining 

S**(x a } c, C) = L(S\ (, = L(x a l a (C 0, c, 0, (20) 

we see immediately that it is a solution of the eikonal equation depending 
on two free parameters. [Note the duplication of notation which arises from 
the different meanings to the same object; S* = ub = x a l a ((, ().} By putting 
the requirement on ( and ( that QS** = dS** = 0, i.e., 

L(S*,(,C)x a m a ((,() + BL(S*,(,0 = 
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L(S*,(,Ox a fh a ((,0+ 9L(£*,C,C) = 0, (21) 



where L = ds*L, we can solve for ( and ( in terms of x a , i.e., equations (pl|) 
give us 

C = T{x a )and( = T(x a ) (22) 
except at the caustics when[|5| 

d 2 S** 5 dS** 
5 dS** d 2 S** 

This issue will be discussed in Sec. 5. 

Finally replacing ((, () in (|20|) and differentiating we find 

d a S**(<r(x a ),r(x a )) = L l a (T(x a ),T(x a )). (24) 
Therefore, the function 

S**{x a ,T{x a ),T{x a )) (25) 

satisfies the eikonal equation, and by construction ( p0[) it is adapted to the 
characteristic data defined by the function L(ub,(X) a ^ 3 + - The null nor- 
mals to level surfaces 

S**{x a ,T{x a ),T{x a ) = const. (26) 

are normal to the cuts L(ub,(X) — const, at J + . Note that the fact that 
S**(x a ) is a new solution of the eikonal equation is a consequence of property 
of the eikonal equation of being homogeneous in d a S. 

4 Relation between the Cauchy and Charac- 
teristic Constructions 

In this section we give the connection between the two methods of con- 
struction. Earlier we showed how to relate the Cauchy data, Sc a uch y (x l ), with 
the arbitrary function H(j3, (, () of Sec. 2, so that any solution of the eikonal 
equation can be cast in the form of equation ([T4]). Therefore, there must be a 
relationship of the characteristic construction to the construction via Cauchy 



0. 



(23) 
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data and hence a relationship between the functions L(x a l a ((, (), (, () and 

We first note that though in both methods there is an arbitrary func- 
tion of three variables, in the characteristic method there appear only two 
parameters ((, () while in the Cauchy method there are the three (/?, (, (). 

We can reduce the three to two by solving Eq. (|15|) 

PjTT 

x a l a ((,0-—((3,(,Q = 0, (27) 

for j3 = /3(x a l a ,( ,C) or changing notation and using ub = x a l a , we have 
(3 = (3(ub, C)0- Now thinking of ( p7| ) as an implicit relation defining either 
u B = U{p, C, C) = ff(A C, C) or P = P( u b, C> 0- Note that if we treat H as 
a function of (%, (, C),i- e -j H = H(/3(ub, (, (), C> C) then 

PjTT 

H^d UB H\ a = —(P,U)P (28) 

or 

d/3 



We replace (3 = (3{ub, (,,(,) into the two conditions, Eqs. fllSD and (|17D , 
obtaining 



P(u B X,0x a m a ((,0-m(P(u B ,(,0,(,0\p = o 
P(uB,C,C)x a m a (C,0-^H(P(u B ,C,0,CO\p = 0, (30) 



which appear similar to Eqs.([H]), namely: 



L{u B ,(,Ox a m a {(,0+QL{u B ,(,0\u B = 
L{u B ,CQx a ma((,d+ dL{u B ,(,0\u B = 0. (31) 

We explicitly write \p and \ UB in the 5 -operators to mean that the angular 
derivatives are taken keeping (3 or ub constant respectively; also L means 
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Remark 5 As we mentioned earlier, Eqs. (WW or fl3l\ ) implicitly define ( = 
T(x a ) and ( = T(x a ) everywhere except at the caustics. They can be ap- 
proached in a limiting fashion. 

Given an arbitrary function F((3, (, () and (3(ub, C> C) there is the follow- 
ing relation between differential operators. 



3F(/U,C)l/3 = mPX,0\u B -(dF/d(3)((3,(,Omu B 
9F(/3,C,C)l/3 = QF(PXX)\u B -(dF/dp)(P,CXW\u B - (32) 

Using these relations to replace the 9 and 8 derivative operators at (3 
constant by operators at ub constant in (|3l] ) we obtain: 



(3x a m a {U)-^H{(3XXO\u B + {dH/d(3){(3XX)mu B = 
f3x a m a ((X)-dH(PX,0\u B + (dH/d/3)(f3,(,Omu B = 0. (33) 

Applying relations (j32|) to the function F = ub = x a l a {CX)i thought of 
as ub = U(/3, (, (), via the following steps; 

3F(/U,C)Ls = 3(z a *«) =x a m a (34) 



QF(P,(X)\u = Mu = (35) 

(dF/d(3)d(3\ UB = (du/d(3)d(3\ UB = p- x W\ nB (36) 
we get the following important equation: 

$x a m a = -Q(3\ UB . (37) 
Finally inserting this relation, with Eq. (|29|) , into (|33|) we obtain 

Px a m a ((,()-m(P,(,0\u B -Hx a m a ((X) = 
Px a m a ((X)-m(P,(X)\u B -Hx a m a ((,() = 0, (38) 

which can be rewritten as 
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d(u B (3 - H) \ UB + {u B P ~ H)x a m a ((, () 
Q(u B (3 - H)\ UB + (u B P- H)x a m a {C,X) 




0. 



(39) 



Comparing Eq.(g^) with (plD we see that they are identical when we set 



L(u B , C, C) = u B p(u B , C, C) - H(P(u B , C, 0, C, 0- 

From equation (|27|) we also have that 



(40) 



m b = <9#/d/?, (41) 
(3 = dL/du B . 

We see that the two data functions L(ub, (, C) an d H(j3, (, () are related 
by the Legendre transformation, Eqs. (|40"D and (^I|) . We have finally arrived 
at a very simple and beautiful relation between the two methods. An essential 
property of the eikonal equation for this relationship is that it is homogeneous 
in d a S. 



5 Parametric Description of the Wave- Fronts 

Using the methods described above we can construct a general solution of 
the eikonal equation either for the Cauchy or the (corresponding) character- 
istic data. In the last section we showed that they are simply related by a 
Legendre transformation in the variables u B and (3. Once we have the solu- 
tion of the eikonal equation we can study the geometry of its wave-fronts. In 
particular we are interested in the description of the singularities developed 
by them, namely, its caustics. 



A key step in both methods consists of expressing the free parame- 
ters, e.g., (/3, C, C)j contained in the formalism as functions of the space-time 
points, x a . In many cases the problem of inverting Eq. ([l5| - |T7|) or Eq . (^l|) in 
order to get either ((,() or (/3, (, () as functions of the space-time coordi- 
nates x a , can be a formidable task and at times impossible. However, it is not 
absolutely necessary, since it is possible to give a parametric description of 
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the null surfaces defined by Eq. (|H] ) or Eq. respectively (In the sequel we 
follow the path given in ref. || for the case of the stationary eikonal equation). 

In the Cauchy case we have three parameters ((, (,/3) in the initial data 
for the eikonal equation. See Eq.(|l4|). We introduce the new parameter r 
together with (£, £,/3) by means of the following equation; 

r = p-mS** = x a (n a - l a ) - /T^StfO^COI/J I 42 ) 
and the previous equations; 



x a i a (CX)-— (/U,C)lc,C = o. ( 43 ) 



dH 

~df3 

The four Eqs.fl43|) and (^) can be solved for the coordinates x a in terms 
of (P,r, (,C), using the orthonormality of the null tetrad: 

a oh na a . . mH. ia m a m a 

x a = — r + n a ) - r —)l a - —m a — m a . 44 

op p p p 



Equation (^) is not very convenient for the analysis of the wave fronts 
because the parameter j3 does not have a simple geometric meaning related 
with the null surfaces. On the other hand, as we know, the level surfaces of 
S** = u = const in Eq. (|i~4|) define the null surfaces in which we are interested. 
Therefore, a sensible parameterization will be the one that replaces the (3 with 
the parameter u defined by 

u = (3x a l a - H(P, C, C) = L(x a l a , C, 0- (45) 

Constant values of u label the characteristic surfaces themselves and are dif- 
ferent than ub = x a l a . By changing the parameter f3 in favor of u we are 
switching to the characteristic description which provides a better framework 
to study the dynamics of the wave fronts. 

Remark 6 Note that r = / 9~ 1 53S'** = /S^SSm defines an affine parameter 
along the null geodesies that rule the characteristic surfaces u = const. 
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Instead of performing the transformation from the "Cauchy parameteriza- 
tion" to the new set (u, r, (, () we take a shortcut, and start directly with the 
characteristic approach. Using the notation of the previous sections for the 
characteristic problem the new parameters are determined by the previous 
equations: 

u = u(x a ) = L(x a l a ,(X), (46) 



L(x a l a ,(,()x a m a ((,()+ 3L{x a l a ,(,0 = 
L(A,(,(>% a ((,C) + gL(A,(,0 = 0, (47) 

and the new one defined by r = L~ 1 <3dS** yielding 



5L„ 3L „ _ L „ h SSL , , 

r = x (n a - l a ) + —x m a + —x m a + —x x m a m b + — — . (48) 

The coordinates x a can be written in terms of the four parameters u,r,( 
and ( as 

x a = u B (l a + n a ) + (r + m + $$)/ a - $m a - $m a , (49) 

where 

SL 

$ = -j-, (50) 

and the function ub = x a l a is written in terms of the parameters u,( and 
( implicitly by u = L(x a l a ,(,C), i-e., x a / a = L _1 (m,(^,^). (L _1 denotes the 
inverse function of L.) 

Treating Eq. (^9]) as a coordinate transformation between the natural co- 
ordinates associated with the solution, i.e., the («,£,£, r), and the standard 
space-time coordinates x a , the transformation breaks down when its Jacobian 
vanishes. This is a three-surface in the space-time; the caustic set associated 
with the solution. 

After a lengthy calculation we find that this occurs when 

J= l?' X ' V ;% =r l -<T d o = 0, (51) 
d(u,r,C,Q 
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where 



a = 5$ + $$. (52) 

This is equivalent to Eq.(3.1) of reference [Q. 

There is a simple geometric interpretation of Eq.(^) and (|52"D; the shear 
function a of the congruence of null geodesies that generate the surfaces 
u = constant, with the affine parameter r, is given by [H] 



r 2 -a°a ' 

Therefore, the vanishing of the Jacobian fl5"T| ) implies that the shear of the 
congruence diverges. We regain the expression defining caustics from ref.|| 
in the stationary case, namely, 

r 2 - 5 2 L 8 2 L = 0, (53) 

since a = <5 2 L. 

The form of the metric tensor in the new coordinates is 



ds z = n a bdx a dx b 

du ,1 + 9$ + $<E> 

= 2^{dr + du( : ) 

L L 

- / 5a- Q + $^°-6a- Q + $ r s „, Qa° + $(x° - $cr + $ r , , 
+ d({ )+d({ )} 

- ^W + ^( 2 )-2(r 2 + ff V)^, (54) 

where P = 1 + This line element, corresponding to shearing non- 
stationary null coordinates, defined by Eq.(^) reduces to the one given in 
ref.Q in the stationary regime. As pointed out in this reference it might 
be of interest to use Eq.(|54|) as a background metric in linearized gravity 
for higher order perturbations in problems where gravitational radiation is 
important. 
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6 The Eikonal Equation in Asymptotically Flat 
Space-Times 

In a straightforward manner all our results can be generalized to the 
case of arbitrary curved space-times, and the proofs of all the relation above 
follow basically the same path. We will assume that there is given a system of 
local coordinates x a in an arbitrary curved space-time and a two parameter 
family (sphere's worth) of solutions of the eikonal equation, i.e. Z(x a , (, Q 
such that 

g ab (x a )d a Z(x a , C, ()d b Z(x a , C, C) = (55) 

such that its (null) gradient sweeps out the light-cone at x a as (£, Q range 
over the sphere. 

Such characteristic function S = Z(x a , (, Q are one of the main variables 
of the Null Surface Formulation of General Relativity; they contain all the 
conformal information of the space-time JF|. In the special case of asymptot- 
ically flat space-times Z(x a , (, () can be interpreted either as the light cone 
cut of 3 + of the point with coordinates x a , or as the past light cone of a point 
at 3 + with coordinates (u, C>C)0- 

We take the complete solution /3Z(x a XX)i an d define, in an analogous 
manner to the flat-space construction, 

u = S**(x a } (3, C, = PZ(x a } C, C) - H({3, C, 0- (56) 

On Eq. (|56|) we impose the conditions, equivalent to (|15|),(|16D and (|17D , 
namely, 

QS** = pdZ(x a XX)-QH(P,CO = (58) 
dS** = pdZ(x a ,(X)-m(P,(X) = 0, (59) 



and solve for (j3, (, () (as earlier, always possible aside from lower dimen- 
sional (caustic) regions which can be approached in a limiting fashion) in 
terms of the x a . 
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When these are resubstituted into Eg. (|56|) , S** then becomes a new so- 
lution of the eikonal equation since 

d a S** = (3d a Z. (60) 

As in the flat case, we can determine the arbitrary function H(/3, (, Q in 
terms of corresponding data given on a Cauchy surface E. Suppose that we 
are given a coordinate system (r, x l ) such that r = tq corresponds to our 
Cauchy surface, together with suitable Cauchy data Sc aU ch y (x l ) on E. A 
needed generalization of the relationship ai = OScauchy/ dx l from Sec. 2 is 

dScauchyjx 1 ) (3dZ{x\ Tp, (, Q _ Q ^ 

dx % dx l 
which is to be considered as three equations for the determination of x % 
in terms of (/3, (, (), i.e., x % = X l (f3XX)- When these are inserted into 
Eq. (|56|) at r = r we obtain 

H(P, (, = PZ(X\P, C, C), To, C, - S C auchy(X l (P, (, ()). 

in analogy to the results of Sec. 2. 

The characteristic formulation from Sec. 3 is even simpler. Starting 
with any function L(ub,(X) defined on J + , we obtain a solution to the 
eikonal equation with the given characteristic data by 

u = C, C) = L{Z{x a , C, C), C, 0, (62) 

where (C,C) are functions of the coordinates x a such that the equivalent to 
Eq.(|2T|) hold, i.e. when 

55** = mZ(x a XA) -QL{Z,£,t) =0 (63) 
85** = VQZ[x\ CO -VL{Z,C, 0=0. (64) 

Again the relationship between H(j3,(X) an d L(u,(,C) is given by the 
Legendre transformation 

L(Z, C, = zp(z, C, - H{f3(Z, C, 0, 0- (65) 

with 
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Z = dH/d(5, 
P = OL/dZ. 



(66) 



7 Conclusion 

We have generalized the results of ref . || concerning solutions of the flat-space 
eikonal equation. We saw two different means of giving data and solving the 
eikonal equation: the Cauchy, and the characteristic formulation. Each one 
leads to different methods. The two methods are beautifully related by a 
Legendre transformation, Eqs. (f4*0"D and (f4T|). Moreover, all our results can be 
generalized to the case of curved space-times. The characteristic formulation 
appears to be better for the study of the dynamics of the wave fronts. By 
means of a suitable parameterization we could describe the caustics in the 
wave-fronts, and find a simple geometric interpretation in terms of the shear 
cr of the null congruence generating the wave fronts. 
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